Degenerate ground state in a mesoscopic YBa2Cu307_2: grain boundary Josephson 

junction 
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We have measured the current-phase relationship /(</?) of symmetric 45° YBa2Cu307_a; grain 
boundary Josephson junctions. Substantial deviations of the Josephson current from conventional 
tunnel-junction behavior have been observed: (i) The critical current exhibits, as a function of 
temperature T, a local minimum at a temperature T*. (ii) At T « T*, the first harmonic of /(</?) 
changes sign, (iii) For T < T* , the second harmonic of I{ifi) is comparable to the first harmonic, 
and (iv) the ground state of the junction becomes degenerate. The results are in good agreement 
with a microscopic model of Josephson junctions between d-wave superconductors. 



The most important phenomenological difference be- 
tween the high-Tc cuprates and conventional supercon- 
ductors regards the orbital symmetry of the supercon- 
ducting order parameter. In the cuprates the pair poten- 
tial changes sign depending on the direction in momen- 
tum space according toBH A(?9) = Aq cos 2(i? — 9), where 
-d is the angle between the wave vector and the (labora- 
tory) X-axis, while 6 is the angle between the Cu-Cu bond 
direction of the superconductor and the a;-axis_ This un- 
conventional d-wave sy|mmetry was predictedB and ex- 
perimentally confirmedura to be directly measurable in 
the Josephson effect between a high-Tc and a conven- 
tional superconductor. Another consequence of the d- 
wave symmetry is that mid-gap states (MGS) with en- 
ergy e = should form on the free surface of a d-wave su- 
perconductor if A{'d) has opposite signs on incident and 
reflected electronic trajectories. □ The MGS density must 
be maximal for (llO)-like surfaces and this prediction has 
in fact been-|Confirmed by STM microscopy on YBCO sin- 
gle crystalg3 which revealed the MGS contribution to the 
YBCO tunneling density of states. The presence of the 
MGS is expected to influence in a spectacular way also 
the Josephson effect in junctions between d-wave super- 
conductors with different crystallographic orientations. 
Yet no clear manifestation of the MGS in the Josephson 
effect in such junctions has been observed so far, which is 
a challenge for the concept of d-wave superconductivity 
in the cuprates. 

Moreove]p||-|due to possible applications in quantum 
computing,ErlZl there is substantial interest in Josephson 
junctions and circuits with a doubly degenerate ground 
state. Such a state was predicted in an asymmet- 
ric 45° junction (6*1 = 0° and 02 = 45°, the angles 
01^2 are defined in Fig. 1), since odd harmonics of the 
Josephson current I{ip) = ^„/„sinrt(p are suppressed 



by symmctry.oEl The current-phase relation observed in 
Ref. [l^ indeed showed a substantial contribution of the 
second harmonic 12- However, there is a finite supercur- 
rent flowing along the interface in the ground state of 
asymmetric 45° junctions.tl Therefore they do not lead 
to completely quiet qubits in the sense of Ref. 0. 




FIG. 1. Schematic picture of the RF SQUID. The YBCO 
thin film occupies the gray area. The inset shows an electron 
microscope image of the narrow grain boundary Josephson 
junction. 

Motivated by the search for both, the MGS in high-Tc 
Josephson junctions and a quiet qubit, we have stud- 
ied symmetric 45° junctions (i.e. junctions with 9i — 
—62 = 22.5°). In this paper we report the first direct 
observation of several effects exclusive to such junctions: 
temperature controlled sign change of the first harmonic 
of the Josephson current, a nonmonotonic temperature 
dependence of the critical current, and the development 
of a doubly degenerate ground state of the system. 
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Let us start with a description of the theoretical pre- 
dictions for Josephson junctions between d-wave super- 
conductors. In symmetric short junctions the Josephson 
current density is conveniently described in terms of the 
Andreev levels in the junction aalii 
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where (p is the superconducting phase difference between 
the banks, $o is the magnetic flux quantum, /c^? is the 
Fermi momentum, d is the average separation of the 
Cu02 planes, /(e) is the Fermi distribution function, and 
En iv, I?) is the energy of the n-th Andreev level for an 
electron incident on the junction at an angle with re- 
spect to the boundary normal. At a given there exist 
only two Andreev levels with energies ±£((^,-(9). 

The nature of the Andreev levels changes with the im- 
pact angle i?: (i) For 22.5° < < 67.5°, MGS are 
formed at = whose energy is split by a finite phase 
difference ip across the contact. In this range of impact 
angles e{p) can be qualitatively described by euGsi'p) = 
A(7r/4) sin{ip/2)y/V{T:/A), where <rVi'&) < 1 is the 
angle-dependent barrier transparency.E3 (ii) For ji?! < 
22.5° and 67.5° < < 90°, no MGS are formed at 
if = and the Andreev levels resemble those formed in 
a Josephson junction between s-wave superconductors. 
In this range of impact angles e{p) can be qualitatively 
described by £conv(</3) = A(0)[1 - V{0) sin^ {tp / 2)]^ . 

When inserted into the equation for j{ip), the two 
sets of Andreev levels yield contributions of opposite 
sign to the Josephson current, I{ip) = /mgs(<p) + 
-^conv('p)- Close to Tc, when T :s> Aq{T), we can 
approximately write /conv oc P(0)(Aq/T) sinyj and 
Imgs -V {n / 4) { Al /T) sin tp. For a sufficiently 

large ratio I?(0)/I?(7r/4) the sign of the first harmonic 
at high temperatures is therefore given by the con- 
ventional contribution. Lowering the temperature to 
V'I?(7r/4)Ao(r) < T < Ao(T), /conv saturates to 
Iconv OC I?(0)Ao sin whereas |/mgs| continues to grow 
according to /mgs oc — I?(7r/4)(A§/T) sincp. As a result, 
near T* ~ Ao2?(7r/4)/2?(0) the first harmonic will change 
sign and therefore the second harmonic will dominate the 
current, leading to a doubly degenerate ground state. 

In a symmetric junction e(ip,-d) = e{(p,—i})i3 and 
therefore the total current along the interface is exactly 
zero at any T and ip, contrasting with a finite total cur- 
rent in an asymmetric junction. A more detailed analysis 
shows that in symmetric junctions time reversal symme- 
try is broken and small currents of equal magnitude and 
opposite sign fiow along the interface in the left and right 
superconducting banks. Thus, although they are much 
closer to being quiet than asymmetric junctions, even the 
symmetric junctions are not completely quiet. 

An experimental confirmation of the above theoret- 
ical predictions requires the use of mesoscopic junc- 
tions. In fact, since the grain-Jpoundaries are faceted 
on a length scale ~ 0.1 /imjld I[pi) of macroscopic 



junctions necessarily represents a nontrivial average over 
junctions with different misorientations. However, stan- 
dard transport measurements of the critical current Ic = 
maX(^{/((^)} are possible only at temperatures smaller 
than the energy of the grain boundary Josephson junc- 
tion, ^ $o/c/27r. Thetefore we have used the modified 
Rifkin-Deaver methodllil'EJ which offers a unique possi- 
bility to study I{p) at temperatures T much higher than 
the junction energy. This is achieved by connecting the 
banks of the junction to form a superconducting ring (or 
an RF SQUID) so that the phase difference across the 
junction is controlled by the (large) phase stiffness of the 
ring. 




FIG. 2. The phase angle a as a function of ipdc measured 
at different temperatures for sample No. 1 (a) and No. 2 (b). 
From top to bottom, the data correspond to (a) T— 30, 20, 
15, 10, 4.2, 1.8 K and (b) T= 35, 30, 25, 20, 15, 11, 10, 5, 1.6 
K. The data are vertically shifted for clarity. 

The YBCO thin films of thickness 100 nm were pre- 
pared by laser deposition on 45° symmetric bicrystal sub- 
strates. The RF SQUIDs were patterned in the shape of a 
square washer 3500 /imx 3500 fim with a hole 50 /xmx50 
/im by electron beam lithography (see Fig. 1.). Here we 
present the data on two samples with critical current den- 
sities jc ~ 2.6 X 10"^ A/cm^ for sample No. 1 and jc ~ 400 
A/cm^ for sample No. 2. The estimated Josephson pen- 
etration depth Aj is much smaller than the width of the 
wide junction, wi = 1725 /im, and larger than the width 
of the narrow junction, Ws — 0.7 fiui and Ws = 0.5 /im for 
samples No. 1 and 2, respectively. Thus the behavior of 
the RF SQUID is dictated by the narrow junction only. 
The submicron bridge was formed at a position between 
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the defects of the substrate which are visiblean Fig. 1. 

In the modified Rifkin-Deaver method jl^'lla the sam- 
ple is inductively coupled to a high-quahty parallel res- 
onance circuit [Q = 155 and 165 for samples No. 1 and 
2, respectively) driven at its resonant frequency wq- The 
angular phase shift a between the driving current and 
the voltage across the circuit is measured by a RF lock- 
in voltmeter as a function of the external magnetic flux 
'^dci which is conveniently measured in the dimension- 
less units tfdc — 27r<I>rfc/*&o- The phase difference across 
the junction ip was calculated from the a{^pdc) data us- 
ing the coupling coefficient between the RF SQUID and 
the tank coil k'^ = 2.6 x 10"'^ and 3.6 x 10^'^ for samples 
No. 1 and 2, respectively. After inverting the (f = ip{^pdc) 
function, I{^p) can be obtained from Pf{^p) — fdciv) ~ fi 
where f {tp) = I{ip)/Ic, P = 27rL/c/$o, and L = 80 pH 
is the inductance of the RF SQUID. The-details of the 
experimental method are given elsewhereJ13 
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FIG. 3. I{ip) for sample No. 1 (a) and No. 2 (b). From top 
to bottom at ip/-K— 0.5, the data correspond to (a) T— 30, 20, 
15, 10, 4.2, 1.8 K and (b) T= 20, 25, 30, 35, 15, 11, 10, 5, 
1.6 K. The data in Fig. 3b corresponding to T = 20, 25, 30, 35 
were multiplied by a factor 4 for clarity. 



the negligible high-T values to values comparable to 
at low T, suggesting a common origin of both phenom- 
ena. Furthermore, Fig. 4 shows that close to T*, there 
is a local minimum of the critical current Ic as a func- 
tion of T, which is associated with the sign change of /i. 
These results are in a qualitative agreement with theo- 
retical predictinns for I{ip) of 45° junctions with ideally 
fiat interfaces .EZIeI 

We can reconstruct the free energy F of the junction 
as a function of ip from F{(p) = ($o/27r) d4>I[4>)- The 
result is shown in Fig. 5. Note that for T < 15 K, the 
free energy minimum of the sample No. 2 moves away 
from iy9 = 0, and the F{(f) curve exhibits two degenerate 
minima at 95 = ±iy9o, as observed previously in Ref. |l^ 
on asymmetric 45° junctions, see Fig. 5. 
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FIG. 4. The critical current Ic (triangles) and the harmonic 
components I\ (squares) and I2 (circles) of the Josephson cur- 
rent as a function of temperature for sample No. 2. The figure 
is obtained by the Fourier analysis of shown in Fig. 3b. 
Inset: Theoretical prediction for the temperature dependence 
of jc, ji, and j2 for a junction with D — 0.3 and p — 0.3. The 
current densities are plotted in units of the Landau critical 
current density, the temperature is in units of Tc- 



The measured a{ipdc) curves shown in Fig. 2 exhibit 
local minima at low T when ipdc is a multiple of 27r. This 
is qualitatively different from what has been observed be- 
fore for 45° symmetric grain boundary Josephson junc- 
tions on saraales with Wg > 1 fiTa where no such minima 
were found.Ej We believe that the difference is caused 
by the existence of the bicrystal boundary defects with 
a typical distance ~ 1 /im (see Fig. I), which cannot be 
avoided for large junctions. 

The Josephson current calculated from the measured 
a(ipdc) data is shown in Fig. 3. Note the anomalous form 
of I{(p) at low temperatures. The anomalies in sample 
No. 2 are much more pronounced than in sample No. I. 
We believe this is a combined effect of smaller junction 
cross-sections and higher junction quality, as evidenced 
by the much smaller values of Ic in sample No. 2. In Fig. 4 
we plot the first two harmonics Ii and I2 for the sample 
No. 2. The most striking result is that for T* « 12 K, h 
changes sign. In the same temperature region where Ii 
starts to exhibit a downturn, the value of I/2I rises from 



In Fig. 4 we compare the experimental data with a 
theoretical treatment based on the quasiclassical Eilen- 
berger equations which was introduced in the s-wave 
case in Ref. 18 and will be described in detail elsewhere. 



Within our approach the junction is described by two 
phenomenological parameters, the junction transparency 
V and the roughness parameter < p < 00. The tem- 
perature dependence of Ii, I2, and Ic is fit well by our 
theory with T> = 0.3 and p — 0.3. The theoretical jc 
is reported in Fig. 4 in units of the Landau critical cur- 
rent density jo — kpAQ/^od, where Aq is the gap pa- 
rameter at the interface. The experimental critical cur- 
rent densit;^_7c for 45° junctions is smaller than jj?""^ 
by a factorEj 10~^ — 10"*. Remarkably, the abso- 
lute value of jc is also well described hy V = 0.3 and 
p ^ |0,3. In fact, in 45° junctions Aq can be estimated 
fromlld Aq ^ IcRn ^ 10^^ — 1 meV, and therefore 
Jo/ic""' ^ Ao/A|j""^ - 10-2 - iQ-i. Together with the 
theoretical result jc/jo ^ 10^^, this explains the experi- 
mental ratio jc/ic""^- 
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FIG. 5. Free energy F{ip) as a function of the phase dif- 
ference across the weak link. The zero of energy has been 
set so that F{0) = 0. (a) Asymmetric 45° grain bound- 
ary (Ref. |l^). Top to bottom curves correspond to T = 30, 
20, 15, 10, and 4.2 K, respectively, (b) Symmetric 45° grain 
boundary (present work). Top to bottom curves correspond 
to T =20, 15, 11, 10, 5, and 1.6 K, respectively. 

Before concluding let us point out that a large I2 can 
be expected also in a real-space scenario, in which it is 
postulated that due to surface roughness and/or twin- 
ning of the superconducting banks, spontaaeaus mag- 
netic flux is generated along the interface.ll3'E3 We are 
not aware of any prediction of a sign change of Ji within 
the real-space scenario. But we can rule it out also in 
a different way, on very general grounds: I2 corresponds 
to a term —F2COs2ip in the junction free energy, where 
F2 = ($o/47r)/2. In the real-space scenario, F2 is the free 
energy gain due to the formation of spontaneous mag- 
netic flux along the junction. Taking I2 ~ 120 nA for 
the sample No. 2 at low T, we estimate i^2/fcs ~ 1-3 
K. However, since the associated pattern of phase differ- 
ence fluctuations along the junction is not controlled by 
the external circuit (unlike the average phase difference 
across the junction), thermal fluctuations would smear 
the I2 component at T > F2/kB, in disagreement with 
experiment. 

In conclusion, we have found that symmetric 45° junc- 
tions exhibit doubly degenerate ground states. Therefore 
they are of potential use in superconducting qubit fab- 
rication. Tho-qubits based on asymmetric 45° junctions 
are not quictfl since thefe,exists spontaneously generated 
flux along the interfacecll which should change sign be- 
tween the two different ground states of the junction. The 
qubits based on symmetric junctions might be close to 



being completely quiet, since their spontaneously gener- 
ated flux is unmeasurably small .til The key technological 
question is how to increase the energy barrier between the 
two degenerate minima at tp = ifo- An interesting pos- 
sibility seems to be to increase the barrier transparency 
by an appropriate doping of the grain boundary.Ej 
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